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I A. AV.Irxt * 

Presents a mathematical jet engine where 
the measured transition functions of the engine thrust can 
be simulated in a quasi-st ationary operation as well as in 
acceleration or deceleration schedules. Because of its sim- 
plicity and high fidelity, It is especially suited for rep- 
resenting jet engines in digital simulation programs. 
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A SIMPLE PROPULSION SYSTEM MODEL FOR THE SIMULATION OF 
NONLINEAR DYNAMIC THRUST RESPONSE. 

„ , * 
Gunther Schanzer and Peter Krauspe, TU Braunschweig 


Summary! A simple mathematical jet engine description is presented 
with which the measured transition functions of the engine thrust 
can he simulated in quasi- stationary operation as well as in 
acceleration or deceleration schedules. Due to its simplicity 
along with high fidelity, it is especially suited for. the 
representation of jet engines in digital simulation programs. 


1. Introduction / 

The representation of measured thrust-time sequences of jet 
engines in digital simulation programs, of whose accuracy particularly 
high demands are made whenever the jet engine thrust is the primary 
control quantity of a regulation process, often presents problems. 

For one, the transition functions (figure 1) can be represented in 
approximation only by linear control paths of very high order 
(2C and more) on account of the large ratios of delay time and 
recovery time. For another, the transition behavior of the thrust 
changes non-linear ly as a function of the engine rpm (figure 2). 

Finally the time behavior, because of internal energy limitations 
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of the allowable air mass flux, is dependent on the magnitude of 
the control quantity jump; here one must differentiate bet ween 
operation on the so-called steady-state operating line (quasi- 
stationary operation) and the so-called acceleration/deceleration 
schedule. The adaptation of a linear control path of high order 
to the named condition is very costly Cl] • The description of the 
engine dynamics in the phase plane C2] here offers the possibility 
for a simple and clear representation of the linear and non-linear 
portion. The model parameters used depend on engine- and flight 
conditions, such as pressure ratios, air mass throughput, inlet 
temperature, Mach number, etc. For reasons of simplicity, 
one must knowingly forego to include the physical relations between 
the engine parameters into the model, as is generally customary 
in the simulation of the entire engine [ 3 ]. The model parameters 
of the presented model are gained basically from thrust measurements 
(ground station test or simulation), which apply for the particular 
flight condition. 

2. Symbols 

F thrust 

F thrust change with time (= dF/dt) 

F-^ idle thrust 

F ■„ maximum thrust 
max 

t time 

T^ time constant 

N 2 rpm of the high pressure compressor 
^2 rpm changes with time (= dN 2 /dt) 
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i Figure 1: Transition behavior of 

^g -i l Vrsl ion different jet engines. 
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Figure 2: Jump responses of s 

jet engine (measurements) for 
a) large and b) small thrust 
changes. 
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Figure 3 : Phase diagram of the 

thrust sequence of a jet engine 



measured 

simulated 


Figure 4: Comparison of measured 
and simulated thrust curves . 







3. 


Bases and physical relations 


The dynamic non-linear engine behavior, such- as is presented, 
e.g., in figure 2, can be approximated (figure 3) with good accuracy 
by straight lines at least for the multi-shaft engines investigated 
JT3D, JT9D, CF6-50C in the phase plane F, F. The desired function 
F(t) can, since an F value is assigned to every value of F, be gained 
from a simple integration of F: 


F(t) = Fo + j P(F)dt. 
0 


( 1 ) 


In figure 4 a plot of the thrust-time curves calculated in this / 
way is plotted along with the measurements on which they are based. 

As the result of the additional statement of the rpm Ng of the high 
pressure compressor, the typical dynamic properties of a jet engine 
power plant become more pronounced (figure 5) • If one subdivides 
the phase curves in figure 3 and the assigned jump response to a 
change in gas throttling in figure 5 into subregions, one can find 
a clear physical interpretation for every individual region: 

(I) Starting region 

Increase of the engine speed Ng from idle with nearly constant, 
relatively small changes in the speed a/ 2 .i. . One must remember 

that for turbine- jet-stream engines there is a strong non-linear 
relation between the thrust and the speed. Thus, one must assign, 

‘e.g., to the idle thrust F T an idle speed of nearly 60$ of N 0 : 

£-1 IT13.X 

thus, the entire thrust spectrum from F T to F lies within the 
remaining 40$.of the allowable speed range. In region I one obtains, 
despite the non-linear behavior of the speed, a nearly constant 
value of F. ' 


4 


(II) Transition region to maximum acceleration N 9 ■ 

u* i max 

In the phase plane the slope F/F of the phase curve is constant. 

(III) Region of constant thrust change F max with almost linear 
dependency on thrust and speed 

The acceleration of a multi-shaft jet engine is limited in 
this range by the pump limitation of the high pressure compressor. 

In multi-shaft engines the safe distance from the pump limit is 
guaranteed by engine controllers. Thus, the dynamic behavior and 
particularly the maximum acceleration capability is determined in 
region III by the combined action of engine and controller. 

(IV) Transition region to the stationary final value 

The time behavior corresponds to the behavior of a linear 
system of first order with the characteristic time constant T y 
(figure 6): 

' f W = y- Nom. « " ■ F <*>1 = y- AF(t) • ( 2 ) 

J v 1 v 

The slope of the phase curve in the phase plane thus is given by 

(3) 

Let the characteristic time constant of region IV be T^* If one 
compares the transition regions II and IV one can find a formal 
agreements both sections were given in the phase plane by a straight 
line slope f/F. Let the characteristic time constant T-, be defined 

i *’■ 

in region II analogous to equation (3): 
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Figure 5* Relations of thrust and speed of a multi-shaft 
jet engine . 

A comparison with figure 5 shows that for the jet engine investigated 
here one can use the approximation 

. t,=t 2 =t„ 

These time constants characterize time behavior of the jet engine 
for small characteristic value jumps ("small time behavior") and 
lie in the range from 0.25 to 1 second. However, one must note 
that the model parameters F , T., , T n and T 0 apply each time only 
for a given ambient- and flight condition. The model parameters 
presented here refer to ground stand tests under normal atmospheric 
conditions. As is known from flight measurements, larger values 
for T^ and Tg are obtained at higher altitudes as the result of 
decreasing air density and the reductions in air mass flow resulting 
therefrom. 

With the aid of the small-time behavior one can also describe 
those steady state changes which are characterized by relatively 
small control quantity jumps, which the system can follow without 
large delays. The characteristic time constant T y of the linear 
system can be obtained by tangent formation directly from the 
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Figure 6s Jump response and phase diagram of a linear system 


of first order. 


measured curve of the transition function of the path, and namely, 
as shown in figure 2, from the jump response to a) large as well 
as b) small jumps: 

T a = T b = T v (:6) 


It is essential that the path can take on practically all possible 

• ■ 

intermediate values of F which lie below the maximum value F„„ v 

max 

of the acceleration line. Thus, the steady-state behavior of the 
engine is obtained in the phase plane by the number of all points 
(F, F) within the phase curve. 


'm + if: 




Figure 7: Block switch diagram of the engine model. 
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Figure 8-: Examples of thrust commands in the phase s plane- 

a) 75 80$ (steady-state) 

h) 40 $ 0 % (limited acceleration) 

f 

3- Model description 

The described relations are shown in the block switch diagram 
in figure 7* According to equation (2). one can determine as a 
function of the difference between desired- and actual value of 
the thrust F the initial value F of a linear system with the time 
constant T . If the thus determined point (F, F) lies in the phase 
plane within the surface bounded by the phase curve (figure 8 a), 
one is dealing with a change of state, and the desired value F 
is equal to the result from equation (2); if the point lies outside 
of the phase boundary (figure 8 b), one must put in for F the 
maximum possible thrust change corresponding to the associated 
value of the ordinate of the phase curve. In figure 7 this action 
is achieved by the function of the limiting device B whose charac- 
teristic line is the phase curve of the jet engine. A subsequent 
integration of F in accordance with equation (1) provides the 
desired value ,F(t). The same is true for the operation of the 
deceleration line or in the region of the steady-state deceleration 
of the jet engine. 
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The good agreement of the measured and simulated thrust 
curves, together with the low programming expenditure , seemed 
to make, the presented jet engine model to he very useful. 
Applications with detailed numerical calculations are contained 
in £4] • 
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